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Abstract. It has been proposed recently the existence of a non-minimal coupling between a 
canonical scalar field (quintessence) and gravity in the framework of teleparallel gravity, motivated 
by similar constructions in the context of General Relativity. The dynamics of the model, known 
as teleparallel dark energy, has been further developed, but no scaling attractor has been found. 
Here we consider a model in which the non-minimal coupling is ruled by a dynamically changing 
Q = ,f,4>/\/J, with f{4>) an arbitrary function of the scalar field (j>. It is shown that in this case 
C*^ the existence of scaling attractors is possible, which means that the universe will eventually enter 

these scaling attractors, regardless of the initial conditions. As a consequence, the cosmological 
coincidence problem could be alleviated without fine-tunings. 

(N 

1 Introduction 

Since the discovery of the universe accelerated expansion, this issue has been one of the most 
active fields in modern cosmology. The simplest explanation for the dark energy — the name 
I given to the agent responsible for such accelerated expansion — is provided by a cosmological 

constant. However, this scenario is plagued by a severe fine tuning problem associated with 
<^ its energy scale [HOE]. This problem can be alleviated by considering a scalar field with 

a dynamically varying equation of state, like for example quintessence dark energy models. 
Scalar fields naturally arise in particle physics (including string theory), and these fields can 
act as candidates for dark energy. So far a wide variety of scalar-field dark energy models have 
been proposed, such as quintessence, phantoms, k-essence, amongst many others [Dll]. 

Usually, dark energy models are based on scalar fields minimally coupled to gravity, as 
CO well as a possible coupling of the field to a background fluid (dark matter). Also, a possible 

. . coupling between the scalar field and the Ricci scalar R is not to be excluded in the context 

of generalized Einstein gravity theories. For example, a quintessence field coupled to gravity 
^ and called "extended quintessence" was proposed in Ref. [5j (see also Refs. IQlTj). K-essence 

models non-minimally coupled to gravity were studied in Refs. [U |9]. 

On the other hand, the Teleparahel Equivalent of General Relativity (TEGR) [TT | [T2] 
is an equivalent formulation of classical gravity, in which instead of using the torsionless Levi- 
Civita connection one uses the curvatureless Weitzenbock one. The dynamical objects are the 
four linearly independent tetrads, not the metric. The advantage of this framework is that 
the torsion tensor is formed solely from products of first derivatives of the tetrad. Thus, apart 
from some conceptual differences, TEGR is completely equivalent and indistinguishable from 
General Relativity (GR) at the level of field equations |lHll2j. Recently the teleparallel gravity 
generalization, the so-called F{T) theory \13\ I14j. has attracted much attention. 

In Refs. |15^ I16j it was proposed a non-minimal coupling between quintessence and gravity 
in the framework of teleparallel gravity, motivated by a similar construction in the context 
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of GR. The theory was called "teleparallel dark energy", and its dynamics was studied later 
in \n\ I18j . However, no scaling attractor was found. Here we generalize the non-minimal 
coupling function (j)'^ — )• /{(p) and we define a = f,(f,/\/y such that a is constant for f{(j)) oc (j>^, 
but it can also be a dynamically changing quantity a{(j)). For dynamically changing a we show 
the existence of scaling attractors with accelerated expansion in agreement with observations. 
As long as the scaling solution is an attractor, for any generic initial conditions the system 
would sooner or later enter the scaling regime in which the dark energy density is comparable 
to the matter energy density with accelerated expansion of the universe, alleviating in this the 
fine tuning problem of dark energy [U H] . 



2 Quintessence field in general relativity 



As is well known [SJ [6l [7] j one can generalize quintessence by including a non-minimal coupling 
between quintessence and gravity. The relevant action reads 



S 



(fx- 



-9 



2 k 
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where R is the Ricci scalar, Sm is the matter action, = SnG and c = 1 (we use the metric 
signature convention (+, — , — , — ) throughout). Also, the parameter ^ is a dimensionless con- 
stant and /{(p) (with units of mass^ and positive) is an arbitrary function of (/) responsible for 
non-minimal coupling. In |191 [201 l2T| [22] adopted f ^ cfp' in the context inflation, other forms 
have been considered in |23[ [23] (for a more general form that includes derivatives see I^SJ ) . 
When including a non- minimal coupling in quintessence, the corresponding Priedmann equa- 
tions are preserved. On the other hand, the pressure density and energy density of quintessence 
are changed to 



where Vd 



2 </> ^ - 2 e ^0 + 6 e (/ + 4 e fl) H\ (2) 

p^ = \^'' + V-QiHf^^^-QiH^f{<P), (3) 
dV/d(j) and /^^ 



df /dep. For find p^f, and p^f) we have used the equation of motion 
'(P + ?>H4>-iRf^^ + V^ = Q, (4) 

in the flat 



calculated from the action ([T]) , and the useful expression R = 6 ( if + 2 
Friedmann- Robertson- Walker (FRW) geometry, see [SI [7]. 



3 Teleparallel equivalent to general relativity 

Following [12], we briefly review the key ingredients of teleparallel gravity. The orthonormal 
tetrad components ha{x^) relate to the metric through 

QiJiu = h"^ h!'^ riab (5) 
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where Latin indices a,b,c, . . . are indices running over 0, 1, 2, 3 for the tangent space of the 
manifold, and Greek indices fj,,iy, p, . . . are the coordinate indices on the manifold, also running 
over 0, 1, 2, 3. In teleparallel gravity, the relevant action is 

5 = dx^hT + Sm (6) 

where h = det(/i"^) = ^J—g■ The torsion scalar T is given by 

T^S^^'^TP^,, (7) 



where 

is the torsion tensor 

is the contortion tensor and 



K^^K' {d,h\-d,h-^), (8) 
K% = (T'^p - T'^^ - T/-") , (9) 



Sr = ^ ( ^'^P + ^'p T% - t\) , (10) 
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is the superpotential. Variation with respect to the tetrad fields yields the teleparallel version 
of the gravitational field equation 

h-' d, [h hx sp'^) - r;, 5/^^ -lh:T = ^ hP e; (ii) 

where @^ stands for the usual energy-momentum tensor. These equations are exactly the 
same as those of GR for every geometry choice. In the case of a flat FRW background metric, 
the tetrad choice have the form 

h\{t) = dmg{l, a{t),a{t),a{t)), (12) 



where o(t) is the scale factor. It is an exact solution of the field equation ( 11 ) and it is seen that 
the corresponding Friedmann equations are identical to the GR ones both at the background 
and perturbation levels [H El [Til [E] • 



4 Interacting teleparallel dark energy 



In what follows we consider a non-minimal coupling between quintessence field and teleparallel 
gravity. The action will simply read 



S 



d X h 



(13) 



Also, we consider a possible interaction between teleparallel dark energy and dark matter as 
was studied in [T7j, since there is no physical argument to exclude the possible interaction 
between them. Varying the action (13) with respect to the tetrad fields yields 

1 



+ 2e/(<^) 



d^cl)d''ci)-v{cp) 



6J'-da(PdPcP = QX. (14) 



3 



Therefore, imposing the flat FRW geometry (12) in (14), we obtain the Friedmann equations 
with 

p^ = l<j>^ + V{cP)-6CH^f{(l)), (15) 

the scalar energy density and 

p^ = ^<P^-V{4>) + 4CHl^4> + 6C (?,H^ + 2H) /(</.), (16) 

the pressure density of field. Here we also use the useful relation T = —6H^, which arises for 
flat FRW geometry. 



In the flat FRW background, the variation of the action (13) with respect to scalar field 
yields the field equation 

4> + 5H^ + V^ + 6(l^H^ = -a. (17) 
where a is the scalar charge corresponds to coupling between dark matter and dark energy, 



and it is defined by the relation SSm/S4> = —ha [TKUISD]. Rewriting (17) in terms of and 
we find the continuity equation for the field 

p^ + 3Hp^{l + u^) = -Q, (18) 

whereas for matter 

Pm + 3H Pm{l+UJm) = Q, (19) 

where a;^ = p^/p^ and uJm = Pm/pm = const > are the equation-of-state parameter of dark 
energy and dark matter, respectively. We define the barotropic index 7 = 1 + such that 
1 < 7 < 2. The parameter Q = (pa is the coupling between dark energy and dark matter. 

5 Dynamical system of interacting teleparallel dark energy 

Of particular importance in the investigation of cosmological scenarios are those solutions 
in which the energy density of the scalar field mimics the background fluid energy density, 
P<t>/ Pm = with C a nonzero constant. Cosmological solutions which satisfy this condition are 
called "scaling solutions" and the cosmological coincidence problem can be alleviated in most 
dark energy models via scaling attractors (it is defined later) [U H] • To study the cosmological 
dynamics of the model, we introduce the followings dimensionless variables: 



"=7f^ ^=7f^' " = ^V7(0), A = - — , a = -^. (20) 
Using these variables we define 

^H_^ -377/^ + 3(2-7) x'^ + AV6aCux 3j 
^~ m~ 2 (2en2 + l) + 2 ' ^ ' 



and the Eqs. (18) and (19) can be rewritten as a dynamical system of ordinary differential 



equations (ODE), namely 



X = —\y^ -VQa^u + x {s-3)-Q, (22) 
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u 



a X 



X' = -Ve (r - 1) 



(23) 

(24) 
(25) 



a 



Ve n 



1 \ ax 



(26) 



In these equations primes denote derivative with respect to the so-called e-folding time 
In a. Also, we define 



Q 



kQ 



n 



f2 

J J. 



VV,4 



(27) 



The fractional energy densities Q. = {k^p) /{SH"^) for the scalar field and background matter 
are given by 

n^ = y^ + x^ -2^u^, n^ = l-n^, (28) 



respectively. The equation of state of the field uj^f, 



reads 



y^ + 4J^aCux + 2{l 



s) C 



u 



(y2+x2)-2^u2 

On the other hand, the effective equation of state We// = {Pm + Pep) I (Pm + 



(29) 



is given by 



eff 



■7 2/-(7-2)x^ + 4W-a^'ux + 2M7 



+ 7 



(30) 



and the accelerated expansion occurs for We// < —1/3. 

Here, we should emphasize that the dynamical system of ordinary differential equations 



(ODE) equations (22)-(26) are not a dynamical autonomous system unless the parameters F 



and n are know. From now on we concentrate on exponential scalar-field potential of the form 
y(0) = Voe-^^'l', such that A is a positive constant, that is, F = 1 (equivalently, we could 
consider potentials satisfying A = —V^p/kV « const, which is valid for arbitrary but nearly flat 
potentials |26| [271 EH] ) • In addition to the fact that exponential potentials can give rise to an 
accelerated expansion (a power-law expansion a{t) cx with p > 1), they possess cosmological 
scaling solutions [H H]. On the other hand, for a general coupling function u = k ^ /(</>), with 
inverse function = f~^{u^/K'^), a{(j)) and n((/)) can be expressed in terms of u (this approach 
is similar to that followed in the case of quintessence in GR with potential beyond exponential 
potential, see [29]). In this case, when a and 11 are functions of u, the four-dimensional 
autonomous system (22)-(26) (A = const) reduces to three-dimensional autonomous system 



(22)-(24) with a = a{u). Combining the equations (24) and (26) we find that 



a{u) = C 



2 n(n)-i 



du 



(31) 
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where C is nonzero constant. For example for 11 = 1/2 then a = const 7^ 0. On the other hand, 
if for instance (2Ii{u) — l)/u(x 1/ {u — Uc) then a{u) oc {u — Uc) — ?• and a'^{u)Il{u) — )• 
when u ^ Uc with Uc a dimensionless constant {a{u) const = 0). Also, in this case /((/>) is 
calculated from (20) and f{(j)) — s- u^/k"^ = const when cj) — t- ±00. Therefore, in general, two 
situations may arise; one where a is a nonzero constant and another where a depends on u with 
an asymptotic behavior a{u) — )• a{uc) = when u ^ Uc and Uc a dimensionless constant. In 
any case we have a three-dimensional autonomous system (22)-(24) (with A = const > 0). Once 
identified the autonomous system we can obtain the fixed points or critical points (xc, Vc, Uc) 
by imposing the conditions x'^ = y'^ = u'^ = 0. From the definition (20), Xc, Vc, ""c should be 
real, with i/c^ and Uc > 0. To study the stability of the critical points, we substitute linear 
perturbations, a; — )• Xc + 6x, y ^ yc + ^V-, and u ^ Uc + 5u about the critical point (xc, yo Uc) 
into the autonomous system (22)-(24) and linearize them. The eigenvalues of the perturbations 
matrix A^, namely, //i, ^2 and determine the conditions of stability of the critical points. 
One generally uses the following classification [H H]: (i) Stable node: /ii < 0, /U2 < and 
/i3 < 0. (ii) Unstable node: /ii > 0, /X2 > and /X3 > 0. (iii) Saddle point: one or two of the 
three eigenvalues are positive and the other negative, (iv) Stable spiral: The determinant of 
the matrix is negative and the real parts of /ii, ^2 and /^s are negative. A critical point is 
an attractor in the cases (i) and (iv), but it is not so in the cases (ii) and (iii). The universe will 
eventually enter these attractor solutions regardless of the initial conditions. Now we are going 
to study the autonomous dynamical system (22)-(24), first for a = const ^ and then for 
dynamically changing a(u), such that a{u) — t- a{uc) = when the system falls into the critical 
point {xc,yc,Uc). Also, the fact that the energy density of dark energy (Q,^) is the same order 
as that of dark matter (fim) in the present universe, suggests that there may be some relation 
or interaction between them. Several different forms of the coupling between dark energy and 
dark matter have been proposed. One possibility usually studied is to consider an interaction 
of the form Q = /3 k pm EO] ■ A second approach is to introduce an interaction of the 

form T pm on the r.h.s. of the continuity equations with the normalization of T in terms of the 
Hubble parameter H, i.e. T/H = r, where r is a dimensionless coupling [4j. Here we consider 
only the first possibility (Q = (3 k pm (p), for both, constant a and dynamically changing a. In 
what follows we shall restrict our analysis to the case of /3 > (/3 = case without coupling), 
but it is straightforward to extend the analysis to the case of negative /?. 



6 Constant a 



6.1 Critical points 



From Eq. (27), for the coupling Q = /3 k /Om it is easy to find that 



Q 



(32) 



with given in the equation 



Substituting ( 32 ) in ( 22 ) we obtain the critical points of 



the autonomous system (22)-(|24[) for constant a 7^ 0. The critical points are presented in the 



table 1 , and these are a generalization of the fixed points that were found in \T7\ [18] . In table 
2 we summarize the stability properties (to be studied below), and conditions for acceleration 
and existence for each point. The parameters q± are defined as 



(33) 
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Table 1: Critical Points for the autonomous system (22)-(24) for constant a 7^ 0. 



Name 




Vc 


Uc 








La 














-1 


7-1 


Lb 








2pi 


<li 


7-1 


7-1 


Lc 








9+ 


< 


7-1 


7-1 


Ld 







V- 


1 


-1 


-1 


A 


2XS 


Le 





X 


v+ 
2A£ 


1 


-1 


-1 



Table 2: Stability properties, and conditions for acceleration and existence of the fixed points 
in table 1. 



Name 


Stability 


Acceleration 


Existence 


La 


Saddle 


No 


/3 = 


Lb 


Saddle 


No 


< 0, a > and /3 > 


Lc 


Saddle 


No 


No 


I.d 


Stable node or spiral, or saddle 


All values 


i > a>Oor^<0, a<0 


Le 


Saddle 


All values 


^ > ^ and a > 



and the parameters v± as 

v± = aC±y^^ (a2^-2A2). (34) 

The point La is a matter-dominated solution (ilm = 1) with equation of state type cosmological 
constant uj^ = —1, that exists for all values of A and (3 = 0. Point Lb is a scaling solution and 
exists solely for ^ < 0, a > and (3 > 0, since in this case the condition < O,^ < 1 is satisfied. 
Accelerated expansion never occurs for these two points because We// > —1/3. The critical 
point I.c is define for ^ > 2(3'^/a'^, a < and (3 > 0, but, it is not a physically meaningful 
solution since fi^ < 0. The points I.d and I.e both correspond to dark-energy-dominated de 
Sitter solutions with 0,^ = 1 and a;^ = Ueff = —1- Point I.d exist for ^ > 2)^ jo? and a > 
or < and a < 0. Point I.e exist for > 2)? jo? and a > 0. We note that the points I.d 
and I.e exist irrespective of the presence of the coupling /3, since = 1 in these cases. 



6.2 Stability 

Substituting the linear perturbations, x — >• Xc -|- y ^ Uc + 5y, and u ^ Uc + 5u into the 
autonomous system (22 )-(24) and linearize them, we find the follows components for the matrix 

3 (2 - 7 ) dQ 

2 

S'jXcVc dQ. 



of linear perturbations A4 

-372/2 + 9 (2-7) xl + 8V6a^ Uc Xc 



Mn 



2(2^^2 + 1) 



dx 



M 



12 



Vex 



Vc 



2 ^ m2 + 1 Qy 



\Xc,yc,Uci 



XcVcUci (35) 
(36) 
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(2en2 + l)^ 2^^^ + 

3 (2 - 7) + 2 V6 a ^ lie \/6 A 



21 



22 



1^ 2en2 + l 2 y 

972/2 + 3 (2-7) x2 + 4^/f)a^ticXc VGXxc-Sj 



M 



23 



2(2en2 + l) 2 
GucVcC (7 yc + (2 - 7) a:?) + 4 \/6 a ^ Xc yc 2 a ^ yc 



(2^^2 + 1^ 



2^^2 + 1 



M 



V6 



a 



31 



A^32 = 0, M33 = 0. 



(37) 

(38) 
(39) 
(40) 
(41) 



From the equation (32) we calculate dQ/dx^ and replace in (35)-(37). The eigenvalues of 



the matrix Jv[ are as follows 
• Point La: 



A'1,2 



/^3 



Point Lb: 



3(2-7) -i±v''+ ""^^fSr^ 



/"1,2 



Point Lc: 



37 
2 ■ 



^3 



37 



(42) 



(43) 



3 (2-7) -l±\ l 



3(7-2)^ 



A'1,2 



Point Ld: 



/"1,2 



3 (-1 ± ^l-|«7eT^C-2A2)) 



/^3 



37 



;W3 = -37- 



(44) 



(45) 



Point Le: 
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3 (-1 ± ^l + |a7eT^2^-2A2)) 



^J■l,2 = ^ , //3 = -37- (46) 

The Points La, Lb, and Lc are saddle points in any case, since ^2 < and fi^ > 0. For 
^ < and a < point Ld is a saddle point, but, for 



a 



and a > it is a stable node. Also, when S, > ^ ^ y^l + then /ii and ^2 are complex 

with real part negative and det(A^) = —9a■y^y^ {a"^ ^ — 2 X'^) < for a > 0. Therefore, in 
this case point Ld is a stable spiral. Finally, point Le is a saddle point, since it exists solely 
for ^ > 2 A^/a^ and a > 0. 



So, for constant a > 0, ^ as in (47) and /? = 0, the universe enters the matter-dominated 
solution La with = — 1 and Wg// = for 7 = 1 (non-relativistic dark matter), but since this 
solution is unstable, finally the universe fall into the attractor Ld (dark-energy-dominated de 
Sitter solution) with = 1 and a;^ = ujeff = — 1- On the other hand, for /3 > 0, ^ < and 
a > the system enters in the unstable scaling solution Lb, but finally the solution Ld is not 



reached since it is a attractor only when a > and ^ > as in (47). Also, for constant a no 
scaling attractor was found, even when we allow the possible interaction between teleparallel 
dark energy and dark matter. 



7 Dynamically changing a 

Now let us consider a general function of non- minimal coupling /(0) such that a depends 
on u and a{u) — )• a{uc) = when {x,y,u) — )• {xc,yc,Uc) i{xc,yc,Uc) is a fixed point of the 
system). The field (p rolls down toward ±00 (x > or x < 0) with f{(j)) — )• u^/k'^ when 
{x,y,u) {xcVcUc). 

7.1 Critical points 



The critical points of the autonomous system (22)-(24) (with a replaced by a{u)) are presented 
in table 3. In table 3 we define ^ = 2 ^ + 1 and C = A + /3 > 0. In table 4 we summarize the 
stability properties, and conditions for acceleration and existence for each point. 

Point II. a corresponds to a scaling solution. The condition < < 1 is satisfied for 
P < (2 - 7) if 



3(2-7)^ L , L \ 3(2-7)^ .... 



or 

\2 



8/?^ \ 3(2-7)^ 
On the other hand, for (3 > y^3/8 (2 — 7) it is required merely that 
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0<«<^^. (50) 
The equation of state for this solution is given by 

'^4' = -r- — w — \ \ + 7-1- (51) 

(^-Ij (2/32 (^_ij_3(7-4)7 + 4(;32-3))+2;32 

For ^ ~ 1 ~ 0), and /3 > 0, we have that « 1. On the other when ^ / 1 and /3 = 0, then 
a;^ 7 - 1 = ujm- 

The accelerated expansion occurs for 

«<-^^^4|^<0. (52, 

Therefore the accelerated expansion no occurs in this solution since necessary ^ is positive to 
ensure < < 1. 

Points II. b and II. c exists for ^ > Both points are scalar-field dominant solutions 

iyi^ = 1), but without accelerated expansion because Wg// = 1- 

Point II. d is also a scaling solution. From the condition < Vt^j, < 1 we have the following 
constraint 

The equation of state of field for this solution is given by 

^4> = — — \ — \ +7-1, (54) 

and for 7 = 1, ^ 1 and \ « (5 {P » I) we have that — 1 (De Sitter solution). The 

universe exhibits an accelerated expansion independently of ^ for 

In the case of non-relativistic dark matter (7 = 1) the accelerated expansion occurs for A < 2 /3. 
Point li.e is a scalar-field dominant solution {Q,^ = 1) that exists for 

0<|<^, (56) 
and gives an accelerated expansion at late times for 

i<^- (57) 
Given that ^ > is necessary for the existence of this solution, then = oj^^f > — 1- 
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Table 3: Critical Points of the autonomous system (22)-(24) for dynamically changing a{u) 
such that a{u) — )• a{uc) = and Uc > 0. We define ^ = 2 ^ lif + 1 and C = A + /?. 



Name 




Vc 


Uc 








^eff 


Il.a 


3(2-7) 





Uc 


3(7-2)^+^ ^ 


Eq. ( 


51 


) ^ + 7-1 


Il.b 


-\fi 





Uc 


1 


1 






1 


II.c 







Uc 


1 


1 






1 


Il.d 




y/3Cl+ 1(2-7)7 


Uc 


A? _|_ 37 , 1 


Eq. ( 


54 


) 


A(7-l)-/3 


2C 






C 


li.e 


AS 
V6 


(' - m 


Uc 


1 


3 


- 1 




A!i_i 

3 ^ 



Table 4: Stability properties, and conditions for acceleration and existence of the fixed points 
in table 3. As in table 3 we define ^ = 2 + 1 and ( = X + (3. 



Name 


Stability 


Acceleration 




Existence 


Il.a 


Saddle 


No 


Eqs. ( 


48), ( 


49 


) and ( 


50) 


lib 


Unstable node or saddle 


No 




^ > 


II.c 


Unstable node or saddle 


No 




^ > 


Il.d 


Stable node or stable spiral or saddle 


\ / 2/< 
^ ^ 3(7-l)+l 


37 
AC 


< t < 37 I C 
- ? - AC ^ A 


li.e 


Stable node or saddle 






< ^ < 


6 

A2 



7.2 Stability 

For dynamically changing a{u) such that a{u) — )• a{uc) = 0, the components of the matrix of 
perturbation M. are written as 



M 



11 



3 (3 (2-7) rc2-7y2) 3(2-7) dQ . 



2C 



dx 



M12 = VQXyc 1 ^Uc,2yc,i 



(58) 
(59) 



M 



13 



Q^UcXc {-fVc- C^-l) xl) 2y/Q^r]cUcxl 



+ 



9Q, 
du 



Uc ^ \Xc,yc,Uci 



M 



21 



3 (2 - 7) Xc yc \/6 A yc 



M 



22 



i 2 ' 

3 ((2-7) xg-37j/g) ^/6Axe-37 
2e 2 



(60) 

(61) 
(62) 



11 



Qiucyc{iyl-{2--i) xl) 2^/6^r]cUcXcyc 
-M23 = ^ 1 



M31 = 0, M32 = 0, M 



33 



(63) 
(64) 



Here we define 77, = 



C — dji \U=Uc- 



Using (32) in (58)-(60) we can calculate the components of M for each critical point and 



the corresponding eigenvalues. 
• Point II. a: 

Hi , 37 



• Point Il.b: 



2-7 2 



PVci I3H 3(2-7) 



2-7' 



2-7 



(65) 



m = M2 = ^, ^3 = -V6/3V^ + 3 (2-7). (66) 

• Point II. c: 



W = /^2 = ^, /X3 = V6/3Ve + 3 (2-7). (67) 

Point Il.d (7 = 1): 



m,2 



3 (A + 2/3) 
4C 



• Point II. e: 



-1 ± 



2/3CC + 3J [CXC 
3 (A + 2/3)2 I 




3?7c 
2C 



(68) 



m = ACC - 37, 



/^2 



/^3 



(69) 



2 ' 2 

Since for the point Il.a we have < | < 3 (2 - 7)^ /2 (Eqs. ([48]), (|49]) and ([50j)), then 
it is always a saddle point. Point Il.b can be unstable node or saddle point. For < ^ < 
3 (2 - 7)2 /2 /32 and 7?c < it is unstable node. For r/c > or | > 3 (2 - 7)^ /2 /^^ it is saddle 
point. Point II. c is a unstable node if < ^ < a/A^ and tJc > 0, whereas for rye < or ^ > e/A^ 
it is a saddle point. 
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To study the properties of stability of the point II. d we shall consider the case where the 
background fluid is non-relativistic dark matter (7 = 1). So, we restrict ourselves to positive 
^ as in constraint (53) with 7 = 1, and A < 2 /3 or equivalently C < 3/3. 



The eigenvalues //i and 112 are real and negatives if 

3 c . 3 



where 



A 



52 + 64/3C(A + 2/3)V(5 



32 /3 (2 A 



> 0, 



(70) 



(71) 



with (5 = (-7 - 6 /3 C + • Then, for r/c < and ^ satisfying ^ the point Il.d is a stable 
node. For r/c > and ^ as in ([70|) it is saddle point. On the other hand for 



(72) 



/ii and fj,2 are complex with real part negative, and the determinant of the matrix Ai that is 
given by 

9vc (2/3C'A|2 + 3C2 (A-2/3) ^-9) 

det{M) = ^ ^, (73) 

4C^e 



is negative if in addition t/c < (also in this case /X3 < 0). Therefore for ^ as in ( 72 ) and ?7c < 



point Il.d becomes stable spiral. The constraint (70) is consistent with the physical condition 



(53) (it can be stronger or weaker). On the other hand, the constraint (|72[) is also consistent 



with (53) provided that A < (^/A. In any case (stable node or stable spiral), point Il.d. 2 is a 



scaling attractor with accelerated expansion (A < 2/3). 
Point II. e is a stable node if 

n P 37 6 

°<^<AC<A^' 

and rjc < 0. Since 1 < 7 < 2 then < y^- On other hand, it is a saddle point for 



37 ^ 6 
AC<^<A^' 



(74) 



(75) 



or rjc > 0. The point II. e is also an interesting solution, it is a scalar-field dominant solution 
{Q^ = 1) and can be attractor (Eq. (74)) with accelerated expansion (Eq. ([57|)) and such that 
^4> = ^eff > — 1- 

So, for dynamically changing a such that a{u) — t- a{uc) = 0, the fixed points Il.d and 
II. e can attract the universe at late times and are equally viable solutions. But since we are 
interested in the scaling attractors, now we are going to concentrate on the solution Il.d. It is 
a scaling solution that has accelerated expansion and at the same time is an attractor (stable 
node or stable spiral). For 7 = 1, the accelerated expansion occurs, independently of ^, for 
A < 2 /3 (Eq. (|55|) ) , and the physical condition < < 1 is ensures for ^ > in accordance 



with (53). Also, for ^ in accordance with (70) and r/c < point Il.d is a stable node, or for 



1^ as in (72) and ?7c < it is a stable spiral. When point Il.d is a stable node, the constraint 



(70) can be stronger or weaker that the physical condition (53). Numerically we find that exist 
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2 4 6 8 10 

N' 



Figure 1: Evolution ol ^Im (dashed), (dotdashed), oj^ (dotted), <j^e// (solid) and a{u) (red 
line) with 7 = 1, A = 2.4 and ^ ~ 3.5 x 10~^, in the case where the coupling to dark matter 
/3 changes rapidly from /3i = to /3 = 4.1 (6 = 0.035 and ui = 0.21). We choose initial 
conditions Xi = 5 X 10~^ Vi = 2.6 X lO"*' and m = 7 x 10~^ and by way of example we 
consider the function a{u) = Uc — u with Uc = 1 and tJc = —1. The universe exits from the 
matter-dominated solution La with constant a, = —1, Wg// = and approaches the scaling 
attractor II. d for dynamically changing a{u) with fi^ ~ 0.70, ~ —0.90 and tOgff ~ —0.63. 
We note that a;^ presents the phantom-divide crossing during cosmological evolution. 



A*(/3) < 2/3 such that for A < X*{/3) then A > (/X and the constraint (70) is weaker that 
([53]). For example, for /3 = 0.6 then A*(/3) ^ 0.82, and if A = 0.2 the constraint (70) implies 
18.75 < I < 22.96. Moreover, the physical condition ([53]) implies 18.75 < | < 22.75. Also, 
in this example, for ^ 19.95 then ujfj, w —1.07 (Eq. (54)) and Q^f, ^ 0.70 (with accelerated 
expansion since A < 2/3). However, this is not a suitable range for ^, since from solar system 
experiments (see for instance [6j), necessary ~ 1 (or ^ ^ 0). Finally, when point II. d is 

0, the physical condition (53) is 



4.1 we have that (72) implies 



a stable spiral, that is, ^ in accordance with (72) and r/c < 
ensures whenever A < (^/A. For example, for A = 2.4 and /3 
0.23 < ^ < 2.90, whereas the constraint (53) implies 0.19 < ^ < 2.90. In Fig. 1 we show the 
case when the system approaches the fixed point II. d in the above example for A = 2.4, /3 = 4.1 
and for instance ^ w 1.00064 (^ ~ 3.5 x 10~^ and Uc = 1). In this figure, by way of example, 
we consider the function a{u) = Uc — u with Uc = 1 and rjc = —1. Also, following [U |31], we 
consider a non-linear coupling 



/3(n) 



W-fSi) tanhC^) +/3 + /3i 



(76) 



that changes between a small /3i to a large /3 in order to that initially the system enters in the 
matter-dominated solution La (constant a) with /3i = and eventually approaches the scaling 
attractor Il.d with w —0.90, Q^jj w 0.70 and accelerated expansion (A < 2/3). 



8 Conclusions 

It has been recently proposed |15l [T6] a non-minimal coupling (j)'^ between quintessence and 
gravity in the framework of teleparallel gravity, motivated by a similar construction in the 
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context of GR. The theory has been called "teleparallel dark energy", and the cosmological 
dynamics of the model was studied later in \n\ I18| . but no scaling attractor was found. 

In this paper we generalize the function of non-minimal coupling </>^ — )• /{(p) and we define 
a = /,0/\/7 such that a is constant for /((/>) oc c/)^, but it can also be a dynamically changing 



quantity For constant a > and as in (47) (/3 = 0), the universe enters the matter- 

dominated solution La = 1) with = —1 and We// = 0, but since this solution is 
unstable, the universe will eventually fall into the attractor I.d (dark-energy-dominated de 
Sitter solution) with Q,^ = 1 and uj,f, = coeff = — 1- Also, for constant a no scaling attractor 
was found, even when we allow a possible interaction between teleparallel dark energy and 
dark matter. 

On other hand, for dynamically changing a{(f)), and when it can be expressed in terms 
of u = K sj f{4>) such that a{u) a{uc) = (we notice that {xc,yc,Uc) is a fixed point 
of the system), we show that interesting solutions exist (see table 3), in particular scaling 
attractors with accelerated expansion. These are non-minimal generalizations of the fixed 
points presented in [H |30l [32] for coupled dark energy in GR. The final attractor is either 
the scalar-field dominated solution II. e or the scaling attractor II. d, both with accelerated 
expansion. The fixed point II. d is a scaling attractor that can be a stable node or stable spiral 
when the background fluid is non-relativistic dark matter (7 = 1). When point II. d is a stable 



spiral, that is, ^ = 2 -|- 1 as in constraint (72), and 



_ da{u) . 

Vc ~ , \u=Uc ^ 'J 5 

au 

it is possible to find a range for ^ in accordance with solar system experiments (^ ~ 1) [6j, and 
compatible with cosmological observations. For example, in Fig. 1 the system approaches point 
II. d with ~ 0.70, Qm ~ 0.30, cj^ w —0.90 and accelerated expansion, in agreement with 
observations. Also, we note that lo^ presents the phantom-divide crossing during cosmological 
evolution, as has been stressed in Refs. |15l [TBI HB] for teleparallel dark energy. Since the 
scaling solution II. d attracts the universe at late times regardless of the initial conditions with 
^Icf, comparable to and accelerated expansion, the cosmological coincidence problem could 
be alleviated without fine-tunings [H H] . 
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